Whispering gallery like modes along pinned vortices 



O 

» 

o 
O 

(N 



5-H 



I 

C 

O 

o 



> 

oo 
On 

in 



X 



Piotr Marecki and Ralf Schiitzhold 

Fakultat fur Physik, Universitat Duisburg-Essen, Lotharstrasse 1, 47057 Duisburg, Germany 

(Dated: October 28, 2011) 

Employing Unruh's analogy to gravity, we study sound propagation in stationary and locally 
irrotational vortex flows where the circulation is wound around a long (rotating) cylinder. Apart 
from the usual scattering solutions, we find anomalous modes which are bound to the vicinity of 
the cylinder and propagate along its axis - similar to whispering gallery modes. For supersonic 
flow velocities (corresponding to an effective ergoregion), these modes can even have zero frequency. 
Thus they should be relevant for the question of stability or instability of this set-up. 

PACS numbers: 43.20.+g, 04.80.-y, 67.25.dt, 43.35.+d 



Introduction The full characterization of sound 
modes propagating within a given flow profile is a ma- 
jor problem in fluid dynamics and often reveals very rich 
physics. Even for stationary flows, which admit a separa- 
tion ansatz where the linear perturbations can be labeled 
by their conserved frequency uj, this problem is highly 
non-trivial: For static systems, depending on what sort 
of scenario is considered, the dynamics of perturbations 
is governed by equations of Schodinger id-tip = Hip or 
d'Alembert <9 t 2 $ = P$ type. In such cases the full char- 
acterization of solutions follows from the spectral analysis 
of the differential operators H or T>. For stationary sys- 
tems, however, the following inherent difficulty appears. 
The equations assume the form d 2 <fr-\-A = where 
the two operators A and B do not commute in general, 
and therefore their spectral content has no direct signifi- 
cance for the problem at hand. As a result, questions like 
the completeness of solutions or the existence of unstable 
modes with $s(ui) < are far more difficult to address. 

The wave equation for sound in a locally irrotational 
and stationary background flow has the form mentioned 
above, 9 t 2 3> + Adt& = B$>. Precisely the same struc- 
ture arises for scalar fields propagating in a stationary 
space-time. Moreover, as discovered by Unruh [1], there 
is an exact analogy between the two cases: Let us con- 
sider a fluid with density g and velocity v, whose pres- 
sure p is a function of g only p = p(g), i.e., the fluid is 
barotropic. The perturbations, i.e., sound waves, can be 
parametrized by a single potential <f> via 8v — V$ and 
Sg = g&/c 2 . Neglecting viscosity, they obey the same 
wave equation as a scalar field in a curved space-time 
described by the effective acoustic metric [1] 



ds 2 = - (c 2 s dT 2 -[dR-v dTf 



(1) 



where T, R are the laboratory coordinates and c s is the 
speed of sound c 2 = dp/dg. This analogy allows us make 
use of many geometrical tools and concepts of general 
relativity in fluid dynamics. 

Vortex flow In the following, we consider a station- 
ary and locally irrotational Vxt) = flow around a long 
cylinder. Aligning the coordinate Z-axis with the sym- 
metry axis of the cylinder, we assume the flow velocity 



to be v = v(R)e v in cylindrical coordinates Z, R, ip. The 
condition Vx» = then implies v = e v K/R, where k 
determines the circulation. 

For normal fluids, such a profile approximates the sta- 
tionary flow around a rotating cylinder. Indeed, if the 
fluid is incompressible "V-v = 0, the above velocity profile 
provides an exact solution of the Navier-Stokes equations 
(similar to a tornado away from the core) . While this so- 
lution gets modified in more realistic models of normal 
fluids, for superfluids (such as 4 He II at low tempera- 
tures), vorticity can only occur in the form of vortices 
with a quantized circulation, which is thus topologically 
stabilized and not so easy to create (or destroy) . The flow 
profile v = e^n/ R then corresponds to a vortex which is 
pinned around a long cylinder (e.g., wire), where k is an 
integer multiple of kq = h/M}i e [2]. 

In fluid dynamics, it is often useful to express the prob- 
lem in terms of dimensionless quantities (such as Rossby 
number) in order to exploit the scaling symmetries. Here, 
we do the same and choose the sound velocity at infinity 
Cs(R f oo) = Cqo as reference scale. Using the circulation 
k we get a length £ = k/coq and a time scale T = k/c 2 ^. 
By rescaling the laboratory coordinates T and R, to get 
dimensionless t and r, we also find that (1) becomes 



ds 2 = 



P 



c 2 - 



1 



dt 2 + 2dt dip - dr 2 



2 dp> 2 - dz 2 



(2) 



where c(r) = c s (r)/coo is the dimensionless sound speed 
which approaches unity at infinity. Similarly, the nor- 
malized density reads p(r) = g{r)/ goo- Typically, both 
decrease a bit when approaching the cylinder where the 
velocity v increases and thus the pressure p drops. 

The coordinates t, z and <p have their standard ranges, 
but r is restricted to r € (r w , +oo), where r w is the re- 
scaled wire radius. Note that the acoustic metric (2) 
possesses an ergoregion goo < for small enough radii r < 
1/c, if such are in the allowed range, i.e., if r w < l/c(r w ). 
Appearance of ergoregion means that the flow velocity v 
exceeds the local speed of sound c s somewhere (e.g., near 
the wire). As we shall see, the presence an ergoregion 
has profound consequences for the sound modes. 
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Geometric acoustics According to Unruh's analogy to 
gravity [1], sound modes in the vortex profile are solu- 
tions of the wave equation of a massless scalar field in 
the metric (2). However, before investigating the full 
wave equation, let us get some insight via the WKB ap- 
proximation - which amounts to studying sound rays. 
They are null geodesies in the space-time (2) which can 
be found in complete analogy to the general relativistic 
Kepler problem. We search for sound rays x a (r) , and find 
that the problem is reduced to quadratures due to exis- 
tence of four independent first integrals. The space-time 
(2) admits three symmetries with the Killing vectors dt, 
d v , and d z . Via the Noether theorem, this implies the 
conservation of the energy E = ft(ft + ip), the angular 
momentum J = Q(r 2 ip — t), and the axial momentum 
P = Hlz, where the dot denotes d/dr. Here the notation 
Q = p{r) I c(r) and / = c 2 — 1 /r 2 is introduced for brevity. 
Together with the null ray condition x a x a = 0, we can 
express all velocities in terms of these first integrals, e.g., 
i = (E- J/r 2 )/{ttc 2 ) and <p = (E + fJ)/{VLc 2 r 2 ). The 
remaining radial equation reads 



n 2 f 2 



p2 _ E 2 J 2 c 2 

? 

c c 



2EJ 



+ ^ = E 2 - V cS (r) ,(3) 



where we have introduced the effective potential V c g(r) 
which also contains the term E 2 (l — 1/c 2 ). 

The sound rays can now be classified by the following 
arguments. For rfoo, the effective potential V e ff(r) van- 
ishes. Hence all scattering solutions must have E 2 > P 2 . 
For r I 0, of the other hand, the effective potential V c s(r) 
diverges V e s(r | 0) | - ooforJ^O. Thus rays are 
strongly attracted by the vortex in its vicinity. There is a 
cut-off in r, however, provided by r w (wire radius) where 
the sound rays are reflected. If, due to r w , the potential 
T4ff(r) is monotonically decreasing for all r > r w , then 
only scattering solutions (E 2 > P 2 ) exist. If, however, 
there are local minima of V e g(r) at finite r G [r w ,oo), 
bound rays oscillating around them will exist. As one 
may easily infer from the structure of Eq. (3), this can 
always be achieved by tuning the angular momentum J. 
Choosing, e.g., J = —E, we see that V e ff(r) is strictly 
negative (assuming c < 1 everywhere; see below), and 
that there will exist rays bouncing off r w indefinitely. 
Note that these are counter-rotating rays, i.e., propagat- 
ing against the vortex flow. 

In the special case of constant c and p (i.e., = 1), 
these qualitative arguments can be made precise. The 
maximum of V e s{r) is at r* = y2J/ (J + IE) where 
V^ s (r*) = 0. Let J, E be fixed and r w , P adjustable. For 
r w > r* we only have scattering rays, otherwise there ex- 
ist also bound rays. Generally, for arbitrarily large wire 
radii r w , one can find values of J and E for which r* is 
imaginary, and therefore bound rays exist. For small J, 
the radius r* goes to zero - i.e., bound states require a 
minimum angular momentum. Finally, for large J, the 
radius r* — > \pl which is outside the ergorcgion at r = 1 . 



Surprisingly, in the special case of constant c and p 
(i.e., il — 1) as above, the problem of finding null or- 
bits, r (<p), is exactly soluble in terms of elliptic functions 
(see [7]), the reason being not-more-than quartic depen- 
dence of V e g(r) on 1/r, as is also the case in the general 
relativistic Kepler problem, for example. We note that 
the J 2 /r 4 -term is crucial for the universal existence of 
bound states as discussed above. As we shall see later, 
this remains correct for the full wave equation. However, 
available treatments [4-6] of the subject of propagation 
of waves in vortex backgrounds introduce assumptions, 
which effectively eliminate this term. While this leads 
to the radial equation of simple Bessel form, it also is 
the reason why the family of bound states reported upon 
here is not to be found in the literature [4-6]. 

Wave acoustics In what follows the full wave equa- 
tion will be considered. To this end, let us first discuss 
the field expressions for E, J, and P. They can be ob- 
tained from the pseudo energy-momentum tensor [12] 



(d a $)(d b $) ~±g ab g cd (d c <P)(d d <i>), (4) 

where $ is the velocity potential Sv = V$ of the sound 
waves Sv and g ab is the acoustic metric (2). For each 
Killing vector field £, we get a conserved field quantity 
S = J dS a T ab £ 6 via integrating over the spatial hyper- 
surface dS a . As usual, invariance under time-translations 
leads to the conserved energy, 



E[$] 



d i r^[^ + c 2 [V<S>Y--[d^] 



(5) 



where Minkowski products are implied in the term [V$] 2 , 
i.e., [V$] 2 = (d r $) 2 + (<9 Z $) 2 + (d v ^) 2 /r 2 . The energy 
functional is positive definite as long as c 2 > v 2 every- 
where, i.e., c 2 > 1/r 2 ^. This is not anymore the case 
if the ergorcgion belongs to the spacetime. The two re- 
maining conserved functionals following from the Killing 
vector fields d z and d v , are the axial momentum 



P[$] 



d z <S>. 



(6) 



and similarly the angular momentum J[$] with d z $> — > 
d v <fr. Furthermore, due to the U(l) gaugc-invariance of 
the complexified wave equation, the Klein-Fock-Gordon 
inner product of two solutions 

< ► 



(*i|$ 2 ) 



d A r d, 



$2 



(7) 



with ( &*9a < I ) 2 = ^\d a ^2 — ®2d a &* is conserved, i.e., in- 
dependent of the Cauchy surface over which it is taken. 

Separation ansatz In view of the symmetries of our 
set-up, we approach the problem of solving the wave 
equation by considering modes specified by the following 
ansatz, which reflects the structure of the Killing vectors 



<f>(t, r, ip, z) = 4>{r) cxp{— iut + imp + ip z z} 



(8) 
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For such modes the conserved quantities are related to 
the inner product via E[G>] = cj($|<1>). P[$] = p 2 ($|$) 
and J[$] = m($|$). Since ($|$) and E[$\ are al- 
ways real, solutions with complex frequencies must have 
($|$) = and E[&] = 0. In the absence of an ergo- 
region, E[$] is positive definite and all frequencies are 
real, i.e., the flow is linearly stable [14]. Furthermore, 
the frequency uj and the pseudo-norm (3>|$) of our modes 
have the same sign in this case. As a result, creation and 
annihilation operators are associated to modes with pos- 
itive and negative frequencies, respectively. In the case 
with an ergoregion, the energy can become negative and 
hence this property is no longer true. This can lead to in- 
teresting and related phenomena such as super-radiance 
[8] and the Klein paradox [9]. Since a given frequency 
lj > can be associated to both, creation and annihi- 
lation operators, one can have a mixing of the two and 
thus phenomena like particle creation. 

By inserting the above separation ansatz (8) into the 
wave equation we reduce it to a single ordinary differen- 
tial equation (in radial direction): 



1 d d 

rp dr ^ dr 



1 - 



m c + 2mw 



H4>= {uj 2 -p 2 z )4> = \4>. (9) 



This main equation of our sound-propagation problem 
has several interesting features [7]. First of all, because 
p — > 1 and c — > 1 at r f oo, the solutions <f)(r) at large 
r are either oscillating, for ui 2 > p 2 , or exponentially 
decaying, for w 2 < p 2 . In complete analogy to the ray 
problem, we call the first type of solutions the scattering 
modes, and the second the bound-state modes. Finding 
these modes is then reduced to an eigenvalue problem 
Hcj) = \<fi with H = T> + V, where T> stands for the "ki- 
netic part" involving r-derivatives and V is the effective 
potential. Note that V e s(r) in the sound-ray problem 
corresponds to V on identification of E and J with w 
and m, repectively. For the scalar product 



{4>l\fa} = I dr p{r)r <j>l(r)^ 2 {r) 



(10) 



the operator H is self-adjoint if the Neumann boundary 
condition at r w is assumed 4>'(r w ) = 0, which would just 
reflect the fact that perturbations can not penetrate the 
wire. Even though neither is H the Hamiltonian of the 
problem, nor is the scalar product (10) distinguished by 
the geometry, both of these elements suffice for the fol- 
lowing general statements. First of all, H has bound 
states, A < 0, only if it one can find test functions ip such 
that {^}\%\tp} < 0. The kinetic part, T>, is always non- 
negative, and therefore bound states can only exist if V 
is sufficiently negative. A lower bound on the eigenvalue 
A can be obtained via the minimum of V(r). Then, for 
general profiles of c(r) and p(r) with the aforementioned 



asymptotics, we get the following statements: 

1) Bound states with u> — can only exist if the velocity 
becomes supersonic somewhere c < 1/r, i.e., if the acous- 
tic spacetime has an ergoregion, for otherwise V(u> = 0) 
is non-negative - consistent with £?[$] = cj($|$). 

2) Bound states with to = can only exist if c s be- 
comes sufficiently smaller than c^ near the wire. The 
mechanism for bound states in this case is just the total 
reflection form the region with larger speed of sound - 
which can also occur in a non-rotating fluid. 

3) Independently of the mechanism of total reflection, 
caused by the w 2 term in V, the other terms in V allow 
for bound states with raw > (i.e., co- rotating) only in 
the presence of an ergoregion, i.e., if the flow becomes 
supersonic v(r) > c s (r) somewhere. 

4) For any radius r w > of the wire there are always 
bound states for some values of m and ui. 

These results can be shown in analogy to the sound- 
ray case. Let us first note that the Bernoulli theorem 
(v 2 /2 + h(g) = const for free stationary flow) implies 
that the enthalpy h(g) drops towards the center, because 
v increases. Thus the pressure p and the density g must 
decrease near the wire. Furthermore, for most fluids the 
Griineisen parameter oc dc s /dg is positive (e.g., for He 
II at T = 0, we have dc s /dg G [2.2,2.9] x c s /g, [10]), 
and therefore the speed of sound is also a monotonically 
decreasing function of r. Thus the term w 2 (l — 1/c 2 ) 
in V is negative. For sufficiently large u, one could get 
bound states for m = via total reflection, see point 
2. However, if the change in c is small and thus the 
required frequencies are too large, the underlying fluid 
dynamic description might not be applicable anymore. 

In contrast, the bound states for m ^ mentioned in 
point 4 can occur for smaller values of u and at larger 
length scales. Let us consider the operator H, where in 
the potential V we may vary the parameters to and ui. 
The remaining freedom of choosing p z can be used to 
adjust A. For an arbitrary test function -0, the kinetic 
part {iplVlip} > of the expectation value of H is inde- 
pendent of m and w, but for counter- rotating modes the 
expectation value of V can be made arbitrarily negative. 
For instance for u = —m the expectation value of V is 
negative for all test functions and scales as m 2 . Thus, 
if m 2 is large enough, we get {?/>|'H|V>} < 0, i.e., bound 
states must exist. 

Case of constant p and c The arguments stated above 
prove the existence of the bound states for large enough 
to, but they do not provide information about how large 
to must be for a given set-up and how to depends on 
p z , for example. To get this information, we numerically 
solve Eq. (9) for the special case of constant p and c. 
For moderate velocities v (sufficiently below the speed of 
sound) , this should be a reasonably good approximation. 

In Fig. 1, we plot examples of the dispersion relations 
ui(p z ) of the bound states. In the right plot, we choose 
r w = 5, in which case that the flow velocity docs not ex- 
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FIG. 1: Dispersion relations for bound states. (Left) with er- 
goregion, r w = 0.3; (solid) co-rotating waves m — 5; (dashed) 
counter-rotating waves m = —5. (Right) without ergoregion, 
r m — 5, rotated and zoomed view for counter-rotating waves, 
rn — —3. The region uj \p\ is filled with scattering states. 

ceed Coo/5 and thus the density should be constant up to 
a few percent. In agreement with statements 2 and 3, we 
only find counter-rotating bound states, here we plotted 
modes with m = —3. In agreement with statement 1, 
there is a gap for such modes, i.e., a minimum frequency, 
of u; m i n ps 1.7. Above this gap, the dispersion relation is 
almost linear and very close to the u = p z line - note the 
different scales on the (p — w)/2 and the (p + u))/2 axes 

- which means that these modes propagate with almost 
the sound speed along the vortex. 

For comparison, we plotted an example with an er- 
goregion r w = 0.3 in Fig. 1 (left). Even though such a 
profile is probably hard to realize experimetally (espe- 
cially with q ps const), one could imagine keeping the 
density nearly constant by applying a suitable external 
potential V such that the stationary Bernoulli equation 
reads v 2 /2 + h(g) + V = const. Furthermore, one could 
rotate the cylinder with the same velocity as the inner- 
most layer of the fluid in order to stabilize the set-up 
better. This profile allows co-rotating bound states (here 
in — 5) which have a negative group velocity in addition 
to the counter-rotating modes (here m = — 5). In this 
case, neither co- nor counter-rotating modes have a gap 

- at discrete values of p z , the frequency and thus the 
energy E = w($|$) vanishes. Such zero-energy modes 
are not uncommon in problems involving perturbations 
of vortices [11] and may indicate an instability. 

Conclusions For a vortex pinned at a wire, we studied 
sound propagation via Unruh's analogy to gravity. On 
general grounds, we predict the existence of bound states 

whispering gallery like modes - based on the geometric 
acoustics approximation as well as the full wave equation. 

It should be possible to verify the characteristics of 
the bound states in experiments. For example, let us 
consider superfluid 4 He with a typical speed of sound 
c s = 2.4 x 10 4 cm/s. For a singly quantized vortex we 
get K = 1.6 x 10 _4 cm 2 /s. This leads to the length 
scale £ = 6 x 10~ 9 cm which is smaller than the van 
der Waals radius of Helium, and to the frequency scale 



c oo/ K — 3-6 x 10 12 Hz above the roton frequency. At these 
scales, fluid dynamics (which is the basis of our sound de- 
scription) is not valid anymore. Thus, for observing these 
bound states, it is probably better to trap many circula- 
tion around the wire, which linearly increases k and thus 
£, while decreasing the frequency of the modes. 

We remark that the family of bound states presented 
here is distinct from the phenomenon of Kelvin waves, 
known from normal- and superfluid dynamics. The vor- 
tex considered in this Letter is pinned and not allowed to 
move or be deformed. Accordingly, the dispersion rela- 
tion (e.g., in Fig. 1) of the modes considered here is much 
stiffcr and more sound-like than that of the soft Kelvin 
waves with U){p z ) ~ K,p 2 z \og{p z )/2 for small p z , see [2]. 
In fact, in the absence of ergoregion the bound states 
discussed here do not even exist for small p z , where the 
above Kelvin formula is intended to be applicable. 
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